Seismic diagnostics of red giants: first comparison with stellar models 
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ABSTRACT 

The clear detection with CoRoT and KEPLER of radial and non-radial solar-like 
oscillations in many red giants paves the way to seismic inferences on the structure of 
such stars. We present an overview of the properties of the adiabatic frequencies and 
frequency separations of radial and non-radial oscillation modes for an extended grid 
of models. We highlight how their detection allows a deeper insight into the internal 
structure and evolutionary state of red giants. In particular, we find that the properties 
of dipole modes constitute a promising seismic diagnostic tool of the evolutionary 
state of red-giant stars. We compare our theoretical predictions with the first 34 days 
of KEPLER data and predict the frequency diagram expected for red giants in the 
COROT-exofield in the galactic center direction. 

Subject headings: Stars: evolution — stars: interiors — stars: oscillations — stars: 
late-type 

1. Introduction 

Red giants are cool stars with an extended convective envelope, which can, as in main se- 
quence solar-like stars, stochastically excite pressure modes of oscillation. Although stochastic 
oscillations were already detected in a few red giants from ground and space observations (e.g. 
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Frandsen et al.ll2002l : |Pe Ridder et al.ll2006l:lBarban et al.ll2007h it has been only thanks to the pho- 
tometric space mission COROT ( Baglin et al. 12002 ) that an unambiguous detection of radial and 
non-radial modes in a large number o f red-giant stars has been achieved (|De Ridder et al.l 12009c 
Hekker et al.ll2009 : Carrier et al. 2010 ). That confirmation has opened the way to the seismic study 
of the structure and evolution of these objects that play a fundamental role i n fields such as stel- 
lar ag e determination and chemical evolution of galaxies. The application of iKjeldsen & Bedding 
(11995b theoretical scaling laws, which relate basic seismic observables (the large frequency separa- 
tion - At/- and the fr equency at maximum power -z/ max ) to the stellar global parameters, allowed 
Mosser et al.1 (I2010al) to estimate the masses of CoRoT and KEPLER red giants (adopting a value 
for the effective temperature). M oreover, the combin ation of these scal i ngs w ith the predictions of 
population synthesis models let lMiglio et al.l (|2009ah and lMiglio et al] (|2009b|) to characterize the 
population of CoRoT and Kepler targets. 

These results, showing that a vast amount of information can be extracted from quite easy- 
to-access seismic observables have deeply changed the perception of the predictive capabilities 
of asteroseismology and have strengthened its interaction with the other fields of astrophysics. 
However, much more information is contained in the oscillation spectra of these large number of 
red giants. In this Letter, we analyze the properties of red giant adiabatic oscillation spectra and 
relate them with their evolutionary state. 



2. Stellar Models 



Stellar models were computed with the code ATON3.1 (|Ventura et al.ll2008l . and references 
therein) for masses between 0.7 and 5.0 M Q and different chemical compositions: He mass frac- 
tion, F=0.25 and 0.278, and metal mass fraction, Z=0.006, 0.01, 0.015, 0.02 and 0.03. The energy 
transport in the convective regions was modeled with the classic mixing length treatment with 
®mlt = 1-6. For a given chemica l composition (Z = 0.02, Y = 0.278), models with «mlt = 1-9 
and FST treatment of convection (|Canuto et al.ll 19961) were also computed. The evolution of these 
models was followed from the pre-main sequence until the exhaustion of He in the core for models 
more massive than 2.3M Q , and until the helium flash for the less massive ones. The core He- 
burning phase for low-mass (0.7 — 2.3M ) stars (Red clump stars) has also been followed starting 
from zero age hor izontal branch model s. Microscopic diffusion was not included but its effects on 
red giant models (|Michaud et al.ll2010|) are not relevant for the present study. 



The evolution rate during ascending RGB, descending RGB and core He-burning phases is 
very different and strongly depends on stellar mass. For low-mass stars the time spent ascending 
the RGB may be comparable with that of core He-burning of more massive stars. As a con- 
sequence, observing stars in both evolutionary stages would be equally likely. Concerning the 
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Fig. 1. — Upper panels: propagation diagrams for two models at almost the same luminosity, 
1.5 M Q in the RGB (left) and 2.5M in central He-burning phase (right). Horizontal dotted lines 
limit the solar-like frequency domain for each model. Solid line is Brunt- Vaisala frequency, and 
short- and long-dashed lines correspond to the Lamb frequency for t = 1 and 2 respectively. Lower 
panels: corresponding plots of inertia as a function of frequency for £=0 (black circles), 1 (grey 
triangles) and 2 (dark-grey squares) modes in the solar-like oscillation domain (dashed vertical 
lines). 

internal structure of these models, it is worth mentioning that for a low-mass model (1.5 M , for 
instance) the density contrast (p c / (p), central to mean density ratio) changes from 10 6 at the bot- 
tom of its RGB to 3 x 10 9 at logL/L Q ~ 2. Models in the core He-burning phase, on the other 
hand, have a much lower value of p c j (p) of the order of 2 x 10 7 ; moreover, due to the high depen- 
dence on temperature of the 3a nuclear reactions, they develop a small convective core. At a given 
luminosity, e.g. that of the Red Clump, the value of p c j (p) for a 1.5M RGB model is more than 
10 times that of a He-burning one. So different structures should imply significant effects on the 
oscillation properties. 
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3. Adiabatic oscillation properties 



Adiabatic oscilla tion frequencies were computed with an Eulerian version of the code LOSC 
(Scuflaire et al. hood) ] or models from the bottom of the RGB until a maximum luminosity (log L/L© 2.2- 
3.2, depending on mass) and as well during the phase of core He-burning. In this paper we deal with 
adiabatic computations a nd do not consider the problem of excitation and damping of solar-lik e 



oscillations in red giants (Dzie mbowski et al]|2001c iHoudek & Goughll2002l : Dupret et al.l l2009). 



We use the scaling laws ( Brown et al.ll 1991c Kjeldsen & Bedding! 1995 ) to derive the frequency do 



main in which solar like oscillati ons are expected, and the va lue of the mode inertia as an estimate 
of the expected mode amplitude ( Christensen-Dalsgaardll2004l) . We search oscillation modes with 
angular degr ee I = 0, 1, 2 and j in the domain of frequencies defined by an interval around z/ max 



(Eq. (10) of iKjeldsen & Bedding 



1995). The width of the solar-like frequency domain is taken to 



be 20% larger than the difference between the acoustic cutoff frequency in the stellar atmosphere 

and is mBX . 

The properties of oscillation modes depend on the behaviour of the Brunt- Vaisala (N) and 
Lamb (Si) frequencies. In red-giant models, N reaches huge values in the central regions and 
therefore the frequency of gravity modes (g-modes) and their number by frequency interval (n g ) 
increase with respect to main sequence models. On the other hand, the low mean density makes the 
frequency of pressure modes (p-modes) to decrease. All that leads to an oscillation spectrum for 
red-giants where, in addition to radial modes, one finds a large number of non-radial modes with 
mixed g-p properties. The dominant character of these non-radial modes depends on the separation 
between gravity and acoustic cavities, and may be estimated from the value of the normalized mode 
inertia (E) (see e.g. IChristensen-Dalsgaardll2004l and references therein). Therefore, some non- 
radial modes may be well trapped in the acoustic cavity and behave as p-modes presenting a mode 
inertia close to that of radial modes, while modes with strong mixed g-p character have larger E 
value. Hereafter, we will use the term p-modes in quotation marks to refer to mixed modes with a 
dominant p-character. 

In Fig. Q] we present, in top panels, the £=1,2 propagation diagrams for a RGB 1.5 M 
model (left) and for a core He-burning (He-B) model of 2.5M (right). In the bottom panels we 
plot the variation of the mode inertia with frequency for radial and non-radial modes (£ = 1, 2). As 
mentioned above, the RGB model is ten times more centrally condensed than the He-B one. The 
huge difference in density between the central region and the convective envelope entails a high 
potential barrier between the acoustic and the gravity cavities reducing the interaction between p 
and g modes. As a consequence, we find for RGB models that £ = 1 modes with Ei = \ « E^ =Q are 
quite regularly spaced in frequency. For He-B ones, the coupling between these cavities is more 
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important and £ = 1 modes are mixed modes with Et = \ > Eg =0 . Nevertheless, Eg = i presents 
still a minimum value for modes between two consecutive radial ones showing a somewhat regular 
pattern. Even if the E value is larger than that corresponding to radial modes we can still consider 
those modes, based on the value of E, as observable "p-modes". For £ = 2 modes, the coupling 
between the g- and p-cavities is smaller than for £ = 1 and hence the trapping more efficient. 
Therefore, independently of the central condensation of the model, a spectrum of regularly spaced 
£ = 2 "p-modes" with E e=2 ~ E t=0 is expected. Finally, note that the turning points for acoustic 
modes (tpg defined as the point where z/ max = Si) are inside the convective envelope for the RGB 
model and in the radiative region for the He-B one. 



In the asymptotic theory for p-modes dVandakurovl 1 19671 : iTassoull Il980l : iGoughl 119860 the 
frequencies of two modes of same degree and consecutive order are separated by a constant value 
(Az/) which is approximately independent of I for low-degree modes, and related to the mean 
density of the star. Of course, the asymptotic theory is no longer valid for mixed modes and in 
regions with rapidly varying physical quantities, nevertheless the modes partially or well trapped in 
the acoustic cavity with a dominant p-character ("p-modes") show such regular pattern. We select 
then, as "p-modes" of degree £, the modes with the minimum inertia between two consecutive 
radial modes and use them to compute the large and small frequency separations and to analyze 
their dependence on the stellar parameters and evolutionary state. 



3.1. Large Separation: Az/ 

The mean value of the large frequency separation ((Az/^) = {u n ^ — z/„_i^) averaged over the 
radial order n) decreases as the star ascends the RGB with a denser and denser core and a more 
and more diffuse envelope, but it is not possible to distinguish, on the basis of the (Az/) value 
alone, among different evolutionary states, ascending RGB, descending RGB, or core He-burning. 
Nevertheless, an indirect information about the evolutionary state is provided not by the average 
value (Az/), but by the deviation of Az/ as a function of frequency with respect to its mean value 
(a(Au e )). Radial and £ = 2 "p-modes", as mentioned above, show a very regular pattern, and the 
mean quadratic deviation of Az/(z/) with respect to its mean value over the solar-like frequency 
domain (<r(Az/^)) is always smaller than 5% for all the evolutionary states and masses considered. 
On the contrary, a(Aui) strongly depends on the evolutionary state and while its value remains 
small for more concentrated models, it may get values as large as 50% for core He-burning ones. 

Fig. [2] shows how the scatter of £ = 1 modes depends on Az/ (that is on (p)) and on the density 
contrast. The scatter of £ = 1 frequencies decreases as the luminosity increases from the bottom of 
the RGB (low p c j (p), or highest (Az/)). As the star goes up the RGB, the £ = 1 modes are better 
trapped in the acoustic cavity and the spectra of dipole modes are more regular. These results are 
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consistent with the observational results obtained by Bedding et al.1 (|2010|) for the first 34 days of 
KEPLER observations: radial and quadruple modes for a sample of 50 low luminosity red giants 
( < 30L Q from scaling laws) show a low scatter, while dipole modes present a significantly larger 
one. The large dispersion shown in Fig.[2]around (Au ) ~ 4 pMz and p c j (p) ~ 2 x 10 7 corresponds 
to models burning He in the core. Following the core expansion, the external convective zone 
recedes, increasing the coupling between gravity and acoustic cavities (see. Fig.Q]) and the mixed 
character of oscillation modes. The spectrum of "trapped" £ = 1 modes is then much less regular. 
It is worth pointing out that, for a given mass (1 M , for instance) at Red-Clump luminosity, while 
RGB and core He-burning models have the same value of (Az/ ), p c / (p) may differ by two orders 
of magnitude. 

The distributions of Au and z/ max for the CoRoT-exo field red giant sam ple show a single domi- 
na nt peak locate d at ~ 4 /xHz and ~ 30/iHz respectively ( Hekker et al.ll2009| ) which was interpreted 
by iMiglio et al.l (|2009ah as being consistent with a population of red giants domin ated by Red 



Clum p stars. According to population synthesis simulations done with TRILEGAL ( Girardi et al. 



20051) for the CoRoT field in the galactic center direction, 70% of the stars at the Red-clump lu- 
minosity are low-mass stars in the core He-burning phase, and therefore, with oscillation spectra 
significantly different from those of the other 30% of stars that, at the same luminosity, are in the 
RGB phase. In Fig. [3] we plot the "ensemble oscillation pattern" that we expect for the CoRoT red 
giant population in the galactic center direction. This diagram was obtained by taking the theoret- 
ical adiabatic spectra corresponding to the stellar models whose stellar parameters are closest to 
those resulting from population synthesis computations. Note the large scatter of £ = 1 modes in 
the region around (Au) ~ 4 pMz, and the regular pattern of £ = and 2 modes. Both characteris- 



tics are in good agreement with the observational results found by LMosser et al.l (I2010b|) . 



3.2. Small separations: 5i/ 02 , Sv 03 and 5v i 



According to the asymptotic theory, the mean small frequency separation ((5^02) = (^no — 
^n-12)) is related to the behaviour of the sound speed (c) mostly in the central regions, and hence 
to the stellar evolutionary state. The representation of (5u 02 ) vs. (Au) is in fact considered as a 
seismic di agnostic diagram allowing to derive stellar mass and age for main-sequence solar-like 
pulsators (|Christensen-Dals gaardll 1 9 8 8|) . The corresponding seismic diagram for red-giant models 
is drawn in Fig.|4]a which shows a linear dependence of (£^02) on (Az/ ) , with a slope that increases 
as the mass decreases. In Fig. [4jb we plot the normalized quantity (5v m) / (Au n ) that in main- 
seque nce stars is known to depend mostly on central physical conditions ([Roxburgh & Vorontsov 
2003J). It is worth noticing that we plotted (5^02) for models with different chemical compositions 
and convection treatment, nevertheless, a predominant dependence on mass and radius appears. 
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For a given mass (8vq 2 ) / (Av ) increases with density contrast i.e. as the star expands with smaller 
values of (v ), decreases as the mass increases, and does not change significantly during the core 
He-burning phase. 

In Fig. HJa and b a vertical dotted line indicates the lower limit of (Av) measured from the 
first 34d of KEPLER mission, an d the dot-dashed line corresponds to the fit (^02) = 0.122(Ai/ ) 
proposed by Beddin g etal.1 (12010 ) for those observations. The comparison between our figures and 
Fig.4 in iBedding et al.1 (|2010l) indicates that the predictions from theoretical models of low-mass 
stars (1-1. 5M ) in the low luminosity part of the ascending RGB are consistent with observational 
data, as the scaling based on z/ max and Av also suggests. In that paper, the authors also report the 
detection of I = 3 modes, and because of the scatter of £ = 1 mode frequencies, they suggest to use 
a new small separation 5v m (n) = 0.5 (von-i ~ 2 v 3n _ 2 + ^nn) instead of the classic <5z/ 13 . Note that 
the observational data provide a ratio 5v 03 /5v 02 ~ 2.2 ( Bedding et alJboiO ) instead of 2 predicted 
by the asymptotic theory. In th e bottom panels of Fig . [4] we plot the theoretical results for 8v 03 and 
its relation with Av and 5v 02 . IBedding et al.1 (|2010h results are represented by a dash-dotted line 
while the dotted line corresponds to our fit Sv 03 ~ 2.42 Sv 02 . The dependence of 5vq 3 on stellar 
mass and radius is similar to that of Sv 02 , therefore no additional information should be expected 
from 8v 03 . 

The small frequency separation 5v i(n) = 0.5 (z/ 0n — 2 v^ n + v 0n+ i) in main-sequence stars 
is also known to be sensitive to the center physical conditions. The asymptotic theory predicts a 
(Su 01 ) = l/3(5v 02 ) relationship. As it is evident in Fig.|4]c the "p-mode" spectrum for red-giant 
models does not follow those predictions. In particular, a large number of models have negative 
or very small values of (5v m) independently of (Sum)- Similar values of (8vqi) have also been 
found in the KEPL ER data (Beddin get al.ll2010|) and in the oscillation spectrum of the C0R0T 
red-giant HR 7349 ( Carrier et al.l 2010 ). We note that the largest concentration of negative/small 
(8u 01 ) values correspond to models ascending or descending the RGB. (5v m ) / {Av) vs. luminosity 
confirms this result. 

Searching for a common characteristic in the structure of these models, we find that while 
ascending and descending the RGB the turning points of t = 1 modes are well inside the convective 
envelope. The steady He-burning models have a shallower convective envelope and the turning 
points of £ = 1 modes are inside the radiative region. Fig.[5Ja shows the variation of (<5^oi)/ (Av) 
with the distance (in acoustic radius r(r') = J Q r dr/c) between the bottom of the convective zone 
(BCZ) and the turning point for a £ = 1 mode with frequency close to z/ max (tpi). In Fig. [5]b we 
highlight the behaviour of (8v 01 ) for low-mass RGB models. The scatter of (Sv i) / (Av) rapidly 
decreases as r ipil — r B cz increases (deep convective envelope) and (5v 01 ) / (Av) takes negative 
values for models in which tpi is well inside the convective envelope. 
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4. Concluding remarks 



We presented the properties of the theoretical spectrum of solar-like oscillations during the 
RGB and core He-burning phases of red giant evolution, and analyzed the behaviour of large and 
small frequency separations computed from modes well trapped in the acoustic cavity of these 
stars. The main results of this global overview are the following: 



Independently of the evolutionary state, 1 = 1 and £ = 3 modes trapped in the acoustic 
cavity have an inertia of the same order as that of the corresponding radial mode and behave 
as "p-modes" with frequencies regularly spaced by (Az/). As a consequence, the scatter of 
£ = 2 modes in the folded echelle diagrams is rather small. 

The trapping of £ = 1 modes in the acoustic cavity depends on the evolutionary state. While 
a regular pattern of dipole modes is expected in more centrally condensed models ascending 
the RGB , the scatter significantly increases for core He-burning ones. Therefore the regu- 
larity of the £ = 1 spectrum could be used to discriminate between different evolutionary 
phases. 

(Suq2) and (5u 03 ) depend almost linearly on the large separation, hence on the mean density 
of the model, with a slope that slightly depends on the stellar mass. 

{Su 01 ) seems to reflect the distance between the £ = 1 turning point and the bottom of 
convective envelope. It takes negative (or small) values if tpi is well inside the convective 
envelope, as it occurs in models ascending or descending the RGB. 

The properties of our adiabatic spectr a are in good agreem ent with those found in KEPLER 
data for low luminosity RGB stars ([Bedding et al.ll2010|) : and also with those for CoRoT- 
exofield red giants. In this study we predict a large scatter of dipole modes for "R ed-Clump" 
stars (Az/ ~ 4uHz. iHekker et al.ll2009l : llVIosser etaDboiOal : llVIiglio et al.ll20Q9ah . as well as 
a regular pattern of £ = and 2 modes for all t he red giants. Both fe atures correspond quite 
well with the properties of the spectra found by lMosser et alj (I2010bf) in the CoRoT-exo field 
red-giants (see Fig.[3]and their Fig. 3). 
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Fig. 2. — Mean squared deviation with respect to the main value of Au e=1 , for RGB and core 
He-B models between O.7M (darkest symbols) and 5.0M Q (lightest symbols), as a function of 
mean large separation for radial modes (a), and as a function of the density contrast (b). Vertical 
dotted and dashed lines correspond to 1M models at log L/L Q ~ 1.68 in the RGB and core He-B 
phases respectively. 
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Fig. 3. — Theoretical oscillation pattern expected for £ = (light grey), 1 (black) and 2 (dark- 
grey) modes from the red giant population in the CoRoT-exofield in the galactic center direction. 
In ordinate the large frequency separation for radial modes, and in abscissa the frequency of (n, I) 
modes are normalized by the large separation. 
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Fig. 4. — Small frequency separations for models with masses between 0.7 and 5 M (each 
grey level corresponds to a different mass, from the darkest for O.7M , to the lightest for 5M ), 
Z=0.006, 0.015, 0.02 and 0.03, F=0.025 and 0.278, a MLT = 1.6 and 1.9. Dotted vertical line indi- 
cates the minimum (Ai/) from the fir st 34-d of KEPLER d ata, and the dot-dashed lines correspond 
to the linear fits for these data from Bedding et al.l (|201ol) . a) {8vq2) as a function of (Az/ ), for 
low-mass models (0.7-2. 3, M Q ) from the bottom of the RGB up to log L/L Q 2.5 — 3, and 2.5-5M 
models in the He-burning phase (Y c between 0.9 and 0.1). b) Same as a) but for the (5u 02 ) / (Az/ ). 
c) (Su 01 ) I (Ai/ ) as a function of ((^02) / (Au ) for the same models as in a) and b). Note the high 
concentration of RGB models with small and negative value of (<5z/ i)/(Az/ ). d) (£^03) / (Ai/q) vs. 
(Au ) for models with Y = 0.278, Z = 0.02, and a M ur = 1.9. e) (5u 03 )/(Au ) vs. (5u 02 ) / (Au ) 
and models used in d). f) (^03) as a function of (^02) and models in d). 
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Fig. 5. — a) (5u i) / (Au ) as a function of the distance in acoustic radius between the bottom of the 
convective envelope r B cz and the £ = 1 turning point r tpl for models with i Z = 0.02, Y = 0.278, 
and «mlt = 1-9 and masses between 0.7 and 5 M . b) zoom for low mass models (0.7-2.3M©) in 
the RGB. Same grey level criteria as in Fig.4. 



